Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  965 1382 (N11E12 L=243)  1852    514        3 ( 867/1199-1117/1596) -36.9 E24-W00 N04-N22
      2  366 1431 (N02E53 L=201)  1821    430        3 ( 236/1359- 461/1475)  14.0 E66-E45 S01-N05
      3 1520 1317 (N23W22 L=277)   353     95        2 (1502/1269-1535/1362) -71.6 W20-W23 N19-N25
      4  244 1102 (S19E54 L=201)  1168    325        6 ( 209/1027- 273/1222) -71.6 E58-E50 S23-S11
      5 1528 1011 (N05W29 L=284)  1457    338        3 (1480/ 828-1593/1157) -40.9 W23-W38 S02-N12
      6  188 1101 (S20E60 L=195)   434    117        2 ( 162/1076- 218/1143) -23.2 E63-E56 S20-S17
      7  672  853 (S25E17 L=238)   839    280        1 ( 542/ 814- 796/ 905)  33.7 E25-E10 S29-S19
      8  276  781 (S38E48 L=207)  1562    423        2 ( 178/ 723- 401/ 817)  21.4 E59-E36 S42-S33
      9 1652  400 (S22W65 L=320)  1989    225        2 (1582/ 345-1700/ 468)  66.0 W59-W69 S26-S17
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y