Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  814 1518 (N21E20 L=304)   689    177        0 ( 766/1435- 844/1596) -48.0 E25-E16 N16-N26
      2 1436 1309 (N17W22 L=346)   825    211        1 (1351/1241-1506/1391) -26.6 W15-W27 N13-N21
      3 1828 1279 (N23W56 L= 20)   950     88        2 (1803/1241-1858/1312) -36.9 W52-W60 N20-N25
      4  139 1045 (S13E68 L=256)  1532    115        2 ( 102/ 989- 174/1090)  45.0 E75-E63 S15-S09
      5  351  957 (S17E45 L=280)   533    192        1 ( 288/ 916- 447/ 980)  15.9 E50-E36 S19-S14
      6 1006  711 (S25W03 L=328)   709    235        0 ( 927/ 636-1049/ 814)  76.0 W00-W05 S30-S18
      7 1717  643 (S17W56 L= 20)   339     39        0 (1698/ 624-1739/ 661)  45.0 W53-W57 S18-S15
      8 1386  615 (S25W30 L=355)   294     94        0 (1336/ 590-1421/ 645)  45.0 W27-W32 S27-S22
      9 1476  414 (S35W47 L= 11)   439     63        0 (1441/ 392-1506/ 444)  63.4 W45-W49 S37-S32
     10  504  321 (S61E49 L=275)   767     57        1 ( 472/ 299- 539/ 343)  11.9 E59-E40 S63-S58
  

    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