Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  889 1861 (N45E45 L=256)   554     59        1 ( 854/1843- 920/1879)  36.9 E47-E41 N42-N47
      2  783 1723 (N31E39 L=262)  1332    163        3 ( 719/1680- 855/1769)  -7.6 E46-E31 N28-N33
      3  464 1591 (N15E54 L=247)  1790    251        5 ( 416/1493- 514/1693) -66.8 E62-E49 N08-N22
      4  825 1545 (N21E26 L=274)   550    109        1 ( 809/1478- 848/1596) -63.4 E28-E24 N16-N24
      5 1660 1251 (N26W35 L=336)  2047    478        0 (1585/1117-1741/1412) -43.0 W26-W42 N19-N34
      6 1810 1222 (N29W50 L=351)   935    114        0 (1765/1161-1843/1255) -26.6 W44-W54 N26-N31
      7  367 1210 (S09E46 L=255)   398     99        0 ( 347/1160- 401/1252)  90.0 E47-E44 S12-S06
      8 1342 1071 (N06W17 L=317)   958    325        1 (1253/ 975-1385/1218)  90.0 W12-W18 N00-N14
      9  223 1086 (S19E57 L=244)   656    150        0 ( 206/1032- 250/1158) -82.9 E59-E53 S22-S14
     10 1871  937 (N14W62 L=  2)  2628    245        5 (1824/ 816-1917/1044) -42.5 W53-W68 N08-N19
     11 1797  977 (N14W52 L=353)   377     42        0 (1778/ 949-1817/1000)  63.4 W50-W53 N12-N16
     12 1771  760 (N02W56 L=357)   925    119        2 (1720/ 731-1818/ 802)  -4.4 W49-W62 N00-N03
     13  765  685 (S30E07 L=294)   333    119        0 ( 718/ 648- 817/ 708)  39.8 E09-E03 S33-S27
     14  321  593 (S49E48 L=253)  2309    238        2 ( 235/ 522- 422/ 670)  -6.8 E60-E34 S52-S44
     15 1414  610 (S17W35 L=336)   277     48        0 (1387/ 588-1440/ 640)  76.0 W33-W36 S18-S14
     16 1506  593 (S15W42 L=343)   511    101        0 (1451/ 577-1565/ 609)  20.6 W38-W46 S16-S13
     17  556  419 (S56E19 L=282)   317     16        0 ( 539/ 395- 574/ 441)  76.0 E21-E17 S57-S53
     18 1261  248 (S43W50 L=351)  1643    163        2 (1181/ 213-1342/ 290)  19.7 W42-W59 S46-S39
  

    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