Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  975 1655 (N29E23 L=338)  1195    400        3 ( 884/1520-1042/1779) -77.2 E28-E19 N17-N39
      2  548 1709 (N21E56 L=305)   632     77        1 ( 521/1678- 579/1743) -71.6 E58-E52 N18-N24
      3 1516 1594 (N44W16 L= 17)  1586    301        3 (1412/1489-1619/1685) -19.3 W05-W26 N39-N47
      4  868 1646 (N25E29 L=332)   719    109        2 ( 813/1617- 917/1681)  39.8 E31-E25 N21-N28
      5 1322 1529 (N31W03 L=  4)   511    132        1 (1288/1464-1340/1582) -54.5 W00-W05 N26-N34
      6 1608 1306 (N27W28 L= 29)   660    272        0 (1517/1199-1674/1392) -21.8 W19-W35 N22-N30
      7  272 1115 (S20E53 L=308)   706    149        2 ( 229/1056- 293/1181) -40.6 E59-E50 S23-S16
      8 1792 1070 (N19W48 L= 49)   601    171        0 (1748/ 998-1827/1140) -24.4 W42-W53 N16-N21
      9  195 1048 (S25E61 L=300)   416     79        2 ( 176/1029- 224/1082) -36.9 E64-E57 S26-S22
     10  417  963 (S26E37 L=324)   363    121        1 ( 362/ 925- 463/ 998)  36.9 E41-E33 S29-S22
     11  964  784 (S22W03 L=  4)   381    140        1 ( 885/ 761-1021/ 810)  39.8 W00-W06 S25-S19
     12 1473  328 (S30W60 L= 61)   551     96        0 (1449/ 289-1501/ 366) -50.2 W57-W63 S33-S27
  

    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