Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  646 1699 (N60E24 L= 24)   502     45        1 ( 621/1674- 672/1719) -26.6 E27-E19 N57-N61
      2  512 1649 (N57E41 L=  6)  2196    228        3 ( 428/1588- 607/1702)   4.4 E54-E28 N52-N61
      3  839 1561 (N44E03 L= 44)   649    124        2 ( 773/1532- 902/1584) -24.0 E09-W01 N41-N46
      4 1285 1535 (N33W29 L= 77)   242     40        1 (1257/1512-1305/1553) -33.7 W27-W31 N31-N34
      5 1048 1525 (N37W12 L= 59)   425    147        0 ( 980/1510-1117/1539)  -5.2 W06-W17 N35-N38
      6 1148 1305 (N21W13 L= 60)   723    206        0 (1092/1251-1226/1337)   0.0 W08-W17 N18-N23
      7  241 1270 (N32E61 L=347)   420     77        0 ( 209/1252- 282/1285)   5.2 E65-E54 N30-N32
      8  747 1095 (N15E16 L= 32)   592    167        1 ( 707/1020- 771/1170) -78.7 E17-E14 N10-N20
      9  270 1048 (N17E54 L=353)  1114    201        1 ( 242/ 948- 312/1119) -55.0 E57-E50 N11-N21
     10 1873  989 (S13W63 L=110)   480     67        1 (1853/ 950-1888/1017) -59.0 W60-W65 S15-S10
     11 1798  845 (S20W53 L=101)  1402    302        1 (1728/ 734-1854/ 971)  40.6 W45-W60 S25-S13
     12 1562  869 (S12W32 L= 79)   364    113        1 (1533/ 829-1594/ 930)  45.0 W29-W34 S14-S09
     13 1557  700 (S23W30 L= 78)  1156    227        0 (1463/ 646-1629/ 773)  15.9 W22-W36 S25-S19
     14 1107  676 (S16E01 L= 47)   662    169        1 (1075/ 600-1153/ 736)  63.4 E03-W02 S20-S12
     15 1559  595 (S31W30 L= 78)   821    260        1 (1456/ 531-1661/ 662)  12.5 W21-W39 S33-S27
     16 1187  524 (S28W02 L= 50)   621    125        0 (1152/ 478-1222/ 571) -69.4 W00-W03 S31-S23
  

    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