Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1604 1570 (N23W58 L= 39)   351     70        1 (1575/1553-1641/1583) -31.0 W55-W61 N20-N24
      2 1344 1286 (N14W26 L=  7)  1879    524        3 (1238/1150-1523/1464) -71.6 W21-W33 N03-N27
      3  692 1287 (N30E15 L=326)  1193    375        3 ( 554/1216- 851/1374) -36.4 E25-E05 N22-N38
      4 1807 1228 (S02W59 L= 40)   443     17        0 (1785/1210-1829/1246)  11.3 W56-W61 S03-N00
      5  952 1126 (N15E02 L=339)   403    252        0 ( 862/1056-1055/1202)  17.1 E08-W05 N12-N16
      6 1822  960 (S18W54 L= 35)   383    150        1 (1786/ 906-1858/1027)  29.1 W49-W59 S20-S14
      7  333  978 (N17E48 L=293)  1622    311        3 ( 196/ 947- 468/1018)  -4.2 E62-E35 N15-N20
      8  634  912 (N09E26 L=315)   893    212        2 ( 561/ 871- 723/ 965)   5.2 E30-E19 N06-N10
      9 1162  569 (S23E02 L=339)   212     73        0 (1125/ 554-1198/ 583) -36.9 E04-W00 S24-S21
     10  373  516 (S11E59 L=282)   403     30        0 ( 360/ 489- 390/ 537)  31.0 E61-E56 S13-S09
     11 1199  453 (S32E03 L=338)   338     61        0 (1176/ 423-1218/ 487)  36.9 E05-E01 S33-S29
  

    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