Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  890 1708 (N30E30 L=140)  1890    401        1 ( 813/1587- 985/1780)  82.4 E35-E22 N21-N36
      2 1306 1661 (N40E01 L=169)   673    156        1 (1243/1600-1345/1714) -29.1 E06-W03 N36-N42
      3  429 1660 (N17E63 L=107)  1597    186        5 ( 358/1630- 515/1693)  14.9 E72-E54 N14-N20
      4  494 1474 (N06E45 L=125)   858    126        1 ( 429/1441- 547/1512)  50.2 E49-E41 N03-N09
      5 1808 1403 (N40W50 L=220)   727    105        1 (1775/1354-1838/1440) -23.2 W45-W54 N37-N43
      6  586 1220 (S06E30 L=139)   244     42        1 ( 557/1201- 607/1241)  45.0 E32-E28 S07-S04
      7 1518 1192 (N15W24 L=194)   255     69        0 (1504/1153-1531/1228)  90.0 W23-W25 N13-N17
      8  889 1125 (S05E10 L=160)   251     95        0 ( 876/1077- 907/1170) -68.2 E11-E09 S07-S02
      9 1704 1006 (N11W42 L=212)   609    187        1 (1663/ 915-1741/1088) -45.0 W37-W47 N06-N15
     10  848  870 (S20E06 L=164)   544    212        0 ( 759/ 794- 898/ 920) -20.0 E12-E01 S23-S18
     11  294  830 (S35E49 L=121)   928    189        0 ( 234/ 777- 336/ 907)  39.8 E56-E44 S39-S29
     12  463  759 (S36E31 L=138)   806    210        0 ( 427/ 676- 524/ 849) -45.0 E34-E25 S40-S30
     13 1598  719 (S08W43 L=213)   533    108        0 (1561/ 682-1621/ 769)  90.0 W41-W45 S10-S04
     14 1466  448 (S26W49 L=218)   851    200        1 (1376/ 418-1547/ 494)  48.4 W42-W52 S29-S20
     15 1111  335 (S45W33 L=202)   404     85        0 (1089/ 298-1142/ 372) -56.3 W30-W35 S47-S41
     16  670  309 (S64E04 L=166)  1448    170        2 ( 595/ 268- 738/ 360)   0.0 E15-W06 S67-S60
     17 1568  332 (S28W69 L=239)  1785    165        3 (1490/ 304-1635/ 359)  11.3 W59-W84 S31-S25
     18 1246  228 (S46W54 L=224)  2926    253        6 (1163/ 185-1351/ 299)   3.9 W39-W69 S50-S42
  

    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