Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  690 1528 (N34E27 L= 11)   748    246        3 ( 604/1466- 759/1586) -26.6 E36-E20 N29-N37
      2 1252 1324 (N22W14 L= 52)  1550    372        4 (1080/1230-1419/1430) -27.6 W02-W25 N16-N28
      3  757 1387 (N24E20 L= 18)   329     57        2 ( 727/1363- 794/1409)  26.6 E21-E17 N22-N25
      4  333 1340 (N19E52 L=346)   699    128        1 ( 315/1282- 354/1398) -60.3 E54-E49 N15-N22
      5  940 1269 (N18E06 L= 32)   614    173        0 ( 880/1192-1010/1342) -45.0 E10-E01 N13-N22
      6  782 1253 (N16E16 L= 21)   541    202        4 ( 745/1157- 817/1340) -74.7 E19-E13 N10-N21
      7  426 1306 (N18E43 L=355)   300     99        0 ( 392/1268- 465/1334)  38.7 E45-E40 N14-N19
      8  535 1136 (N08E32 L=  6)   830    239        3 ( 495/1041- 594/1232) -54.0 E36-E27 N02-N13
      9 1223  257 (S51W24 L= 62)   562    139        1 (1174/ 221-1256/ 284) -26.6 W17-W29 S55-S48
  

    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