Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1525 1520 (N18W47 L= 54)  1255    274        1 (1428/1490-1634/1543) -25.0 W39-W54 N13-N21
      2 1813 1410 (N02W67 L= 75)   959    251        0 (1774/1375-1859/1446)  -3.2 W61-W79 N01-N03
      3 1519 1368 (N10W39 L= 47)   267    109        1 (1466/1356-1561/1388) -15.9 W35-W42 N09-N12
      4 1605 1169 (S03W38 L= 45)   753    171        1 (1574/1091-1656/1256) -90.0 W36-W40 S08-N02
      5 1890 1057 (S19W60 L= 67)   529    144        2 (1870/1019-1921/1093)  20.6 W56-W64 S20-S17
      6  159  968 (N23E63 L=305)   867    140        1 ( 122/ 940- 185/ 990)   7.1 E67-E59 N20-N23
      7 1532  745 (S24W23 L= 30)   243     83        0 (1522/ 704-1543/ 786)  80.5 W21-W23 S27-S21
      8 1694  721 (S32W35 L= 42)   285     51        0 (1674/ 686-1717/ 752)  36.9 W32-W37 S33-S30
      9  199  719 (N07E68 L=300)  1201    185        1 ( 156/ 682- 238/ 751) -32.0 E72-E63 N04-N09
     10  404  458 (S12E59 L=308)   758    119        1 ( 377/ 401- 428/ 508)  74.1 E61-E57 S15-S08
  

    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