Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  928 1649 (N47W14 L=161)   355    134        0 ( 892/1597- 972/1693)  20.0 W08-W19 N44-N49
      2 1023 1392 (N28W10 L=157)   465    156        0 ( 946/1337-1070/1430)  16.7 W03-W13 N26-N29
      3 1157 1243 (N16W13 L=160)   192    106        0 (1105/1226-1199/1269)  -9.5 W10-W16 N15-N17
      4 1889 1169 (S09W66 L=213)  2693    217       12 (1825/1105-1936/1257)  17.4 W56-W74 S13-S04
      5  705 1196 (N24E15 L=132)   763    210        1 ( 617/1140- 797/1235)   4.1 E22-E08 N22-N25
      6  413 1104 (N25E38 L=109)   508    148        1 ( 355/1065- 467/1142)   5.7 E43-E33 N23-N26
      7 1502  836 (S15W24 L=171)  2031    607        4 (1307/ 710-1642/1063)  20.8 W08-W38 S19-S06
      8 1831  946 (S20W53 L=200)   632    147        0 (1797/ 887-1856/1021)  36.9 W48-W56 S22-S15
      9  529  607 (S06E41 L=107)   763    239        0 ( 430/ 552- 625/ 679) -45.0 E46-E35 S11-N00
     10 1549  612 (S31W24 L=171)   597    181        0 (1481/ 573-1612/ 675)   9.5 W17-W29 S32-S28
     11 1142  512 (S26E06 L=141)   774    244        2 (1066/ 429-1213/ 562)  17.1 E13-W00 S29-S23
     12  814  417 (S24E30 L=118)   577    188        1 ( 780/ 343- 855/ 493)  36.9 E35-E23 S28-S19
     13 1135  424 (S32E09 L=138)   252    103        0 (1099/ 395-1155/ 471)  33.7 E12-E06 S33-S28
     14 1162  387 (S35E09 L=139)   886    198        3 (1088/ 324-1222/ 470)  33.7 E15-E02 S38-S30
     15  657  368 (S24E43 L=104)   849    202        1 ( 610/ 300- 718/ 453)  84.8 E46-E38 S29-S18
     16 1427  373 (S46W10 L=158)   658    137        1 (1378/ 334-1466/ 428)  40.6 W05-W14 S48-S41
  

    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