Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1595 1696 (N25W70 L= 40)  1015     67        3 (1560/1665-1622/1721)  59.0 W65-W76 N22-N27
      2 1378 1581 (N27W41 L= 11)   591    192        0 (1279/1555-1463/1601) -36.0 W35-W46 N22-N30
      3 1323 1475 (N23W32 L=  2)   732    213        1 (1218/1430-1407/1507) -35.0 W25-W35 N17-N26
      4 1055 1356 (N24W11 L=341)   654    256        3 ( 999/1281-1123/1470) -90.0 W09-W12 N18-N32
      5  295 1427 (N50E53 L=276)   628     64        1 ( 271/1397- 324/1454)  45.0 E58-E49 N47-N51
      6  518 1302 (N36E26 L=304)  1049    203        4 ( 437/1256- 614/1353) -25.0 E33-E18 N31-N40
      7  708 1197 (N24E14 L=315)   308     92        0 ( 669/1165- 750/1233) -63.4 E16-E12 N21-N27
      8 1216 1193 (N11W15 L=345)   238     77        0 (1180/1159-1245/1229) -68.2 W14-W16 N08-N13
      9 1602  969 (S13W33 L=  3)  1051    291        1 (1541/ 853-1656/1104)  40.2 W26-W39 S17-S06
     10  830 1013 (N11E11 L=318)   604    199        0 ( 771/ 945- 901/1092) -70.0 E13-E08 N05-N16
     11  291  946 (N19E51 L=278)  1230    261        3 ( 229/ 859- 329/1067)  50.7 E59-E46 N14-N25
     12  536  946 (N14E31 L=299)   346     82        2 ( 500/ 918- 567/ 976) -59.0 E33-E29 N11-N16
     13 1613  700 (S30W29 L=358)  1018    142        2 (1561/ 648-1676/ 760)  24.0 W24-W34 S32-S26
     14 1409  713 (S22W14 L=344)   314    159        0 (1339/ 688-1468/ 756)   0.0 W09-W18 S22-S20
     15  248  724 (N07E63 L=267)   457     63        0 ( 210/ 710- 281/ 736)  -6.3 E67-E58 N05-N07
     16 1164  580 (S23E04 L=326)   301    138        0 (1118/ 531-1187/ 634)  39.8 E08-E01 S24-S19
     17 1066  305 (S39E22 L=308)   652    139        0 (1038/ 245-1087/ 370)  48.0 E27-E18 S43-S33
     18  825  329 (S30E36 L=294)   588     94        2 ( 789/ 296- 863/ 369)  12.5 E39-E30 S31-S27
  

    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