Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1667 (N33W43 L=197)   810    268        1 (1220/1634-1447/1700) -19.4 W34-W51 N29-N35
      2 1640 1633 (N20W64 L=218)  2311    389        4 (1541/1592-1728/1667) -26.6 W56-W75 N15-N25
      3 1223 1531 (N29W28 L=182)   455    108        0 (1193/1486-1260/1578)  81.9 W26-W29 N25-N32
      4 1660 1479 (N12W55 L=210)   364     59        0 (1634/1449-1688/1512) -78.7 W54-W56 N09-N14
      5 1622 1427 (N11W50 L=204)   491    154        0 (1606/1352-1641/1492)  58.0 W47-W52 N06-N14
      6  702 1284 (N30E12 L=142)   556    148        2 ( 637/1255- 766/1321) -32.0 E16-E08 N26-N33
      7  796 1201 (N22E08 L=146)   343     67        0 ( 769/1168- 829/1242) -80.5 E09-E06 N19-N25
      8 1898 1127 (S15W64 L=218)  1637    236        1 (1855/1050-1941/1188) -70.0 W59-W67 S21-S10
      9 1814 1107 (S13W53 L=207)   341     83        2 (1793/1072-1831/1130)  21.8 W50-W55 S14-S11
     10  224  926 (N19E57 L= 97)   403     79        0 ( 215/ 885- 237/ 972)  71.6 E58-E55 N16-N22
     11 1803  864 (S26W47 L=201)  1692    255        1 (1773/ 739-1833/ 966)  70.3 W43-W51 S35-S20
     12  191  814 (N13E65 L= 89)  1944    380        3 ( 111/ 770- 235/ 875) -10.0 E76-E58 N10-N16
     13  951  826 (S02E09 L=145)   249     78        1 ( 939/ 781- 971/ 867)  63.4 E10-E07 S04-N00
     14 1239  723 (S16W04 L=158)   923    266        0 (1137/ 593-1308/ 846)  33.7 E05-W10 S21-S10
     15  733  680 (S05E26 L=128)   992    215        3 ( 688/ 592- 791/ 776)  79.7 E27-E23 S11-N00
     16 1433  597 (S30W13 L=167)   884    342        4 (1342/ 513-1529/ 745)  29.4 W05-W22 S33-S23
     17 1147  401 (S34E12 L=142)  1353    256        4 (1066/ 320-1241/ 492)  21.3 E20-E02 S38-S29
     18  555  330 (S23E54 L=100)   611    142        2 ( 534/ 291- 575/ 365)  32.0 E58-E50 S25-S20
  

    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