Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1493 1728 (N48W56 L=183)   499     50        0 (1473/1701-1513/1754)  20.0 W51-W62 N45-N49
      2  899 1427 (N31E06 L=120)   753    243        2 ( 832/1343- 943/1520) -59.7 E11-E03 N25-N37
      3  382 1371 (N28E48 L= 79)   455    132        1 ( 351/1324- 407/1414)  68.2 E51-E45 N25-N31
      4 1276 1330 (N22W19 L=145)   897    273        2 (1202/1251-1360/1371)   5.2 W13-W24 N17-N24
      5  258 1298 (N23E60 L= 67)  1077    133        2 ( 237/1228- 288/1352) -81.9 E62-E56 N19-N26
      6 1831 1273 (N13W65 L=192)   928    101        1 (1800/1240-1866/1318) -63.4 W60-W69 N10-N16
      7  189 1232 (N19E67 L= 60)  1272    102        4 ( 160/1180- 216/1282)  60.3 E71-E63 N15-N22
      8 1573 1203 (N12W38 L=165)   422     56        2 (1547/1170-1608/1228) -36.9 W35-W40 N09-N13
      9 1692 1177 (N09W47 L=174)   281     74        0 (1669/1154-1709/1210) -90.0 W44-W48 N07-N11
     10 1853  947 (S07W62 L=189)   665     67        1 (1836/ 914-1872/ 987) -63.4 W59-W64 S08-S04
     11 1532  837 (S10W32 L=159)   372     74        2 (1513/ 794-1549/ 878)  78.7 W30-W33 S12-S07
     12 1591  757 (S16W37 L=164)  1412    303        6 (1497/ 698-1666/ 842)  28.3 W29-W42 S19-S10
     13  969  713 (S14E05 L=122)   685    263        0 ( 920/ 610-1052/ 813)  56.3 E08-W00 S21-S08
     14 1809  687 (S23W61 L=188)  1555    230        3 (1740/ 602-1866/ 749)  34.7 W54-W68 S28-S19
     15 1532  643 (S23W33 L=160)   699    182        1 (1489/ 592-1576/ 699)  74.1 W29-W36 S26-S19
     16  573  486 (S29E37 L= 90)  2440    414        4 ( 420/ 410- 665/ 606)  21.4 E54-E26 S34-S20
     17  781  507 (S28E20 L=107)   314     80        1 ( 757/ 469- 805/ 538) -78.7 E21-E18 S30-S25
  

    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