Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  901 1851 (N57E34 L= 85)   623     95        2 ( 857/1836- 952/1865)  11.3 E39-E28 N55-N58
      2  809 1693 (N41E29 L= 90)   337     67        0 ( 775/1665- 841/1711) -15.9 E32-E25 N38-N42
      3  693 1626 (N34E35 L= 84)   384     91        0 ( 663/1584- 726/1665) -29.1 E39-E30 N31-N36
      4  976 1627 (N38E13 L=106)   571    142        1 ( 929/1596-1051/1646)  12.5 E16-E07 N35-N39
      5 1459 1400 (N29W25 L=144)  2239    478        0 (1250/1313-1649/1464)  -5.5 W09-W40 N26-N31
      6  433 1107 (S02E39 L= 80)   485    153        0 ( 404/1046- 453/1182)  82.9 E40-E37 S06-N02
      7 1659 1139 (N15W41 L=160)   341    114        0 (1642/1093-1700/1163)  53.1 W39-W44 N12-N16
      8  207 1055 (S08E60 L= 59)   869    217        0 ( 183/ 954- 242/1149)  76.0 E62-E56 S14-S02
      9  665 1065 (S02E22 L= 97)   310     80        0 ( 628/1038- 701/1094)  36.9 E24-E20 S03-N00
     10 1000  765 (S15W02 L=121)   194     99        1 ( 972/ 726-1015/ 805) -76.0 W00-W02 S17-S13
     11 1713  696 (S11W53 L=172)   646    147        0 (1663/ 641-1756/ 760)  63.4 W50-W56 S14-S06
     12 1674  493 (S24W60 L=179)  1612    233        0 (1551/ 469-1760/ 521)  15.9 W48-W72 S27-S21
  

    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