Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1312 1781 (N51E08 L=238)   669     78        1 (1268/1757-1354/1804)  14.0 E11-E03 N49-N53
      2  411 1626 (N13E62 L=184)   966     94        1 ( 377/1591- 449/1659)  68.2 E65-E59 N10-N15
      3  286 1428 (S01E62 L=184)   533     84        2 ( 255/1386- 308/1456) -14.0 E66-E57 S03-N00
      4 1684 1433 (N39W33 L=279)   435     55        0 (1656/1406-1706/1456) -21.8 W30-W35 N37-N40
      5  234 1323 (S08E63 L=183)   486     81        0 ( 217/1287- 248/1369) -59.0 E65-E60 S10-S05
      6  314 1228 (S12E51 L=195)   426    109        0 ( 268/1192- 352/1250)  -6.3 E55-E46 S13-S10
      7  662 1158 (S08E24 L=222)   565    221        1 ( 629/1074- 745/1222)  73.3 E25-E20 S13-S03
      8 1352 1174 (N11W14 L=260)   481    112        1 (1297/1130-1393/1216)  -8.1 W10-W17 N08-N12
      9  149 1165 (S18E69 L=176)  1582    120        1 (  96/1127- 212/1191)  11.3 E80-E60 S20-S16
     10 1010 1077 (S04E02 L=243)   210     59        1 ( 985/1052-1030/1104)  90.0 E03-W00 S05-S02
     11  467  946 (S25E32 L=213)   727    210        2 ( 420/ 861- 546/1018)  71.6 E35-E27 S30-S18
     12 1285  773 (S12W22 L=267)   258    113        0 (1234/ 752-1339/ 789)  45.0 W19-W24 S14-S09
     13 1402  654 (S15W33 L=279)   304     85        0 (1364/ 629-1423/ 690)  78.7 W31-W34 S17-S12
     14  540  537 (S49E20 L=225)  1158    265        2 ( 422/ 508- 651/ 598)   4.8 E33-E09 S51-S44
     15 1668  558 (S11W57 L=302)   398     78        1 (1650/ 514-1692/ 588) -76.0 W53-W59 S13-S09
     16 1232  515 (S28W29 L=275)   472    123        0 (1177/ 495-1295/ 542)  20.6 W24-W32 S29-S25
     17 1539  406 (S23W57 L=303)  1419    217        4 (1442/ 359-1639/ 444)  28.1 W49-W65 S27-S18
     18 1098  263 (S48W39 L=285)   932    170        1 (1026/ 218-1167/ 311)   0.0 W30-W49 S49-S46
  

    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