Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1501 1770 (N35W58 L=172)  1399    332        2 (1447/1727-1550/1817)   7.6 W50-W65 N33-N37
      2  435 1414 (N38E33 L= 82)   292     50        1 ( 411/1387- 464/1433)  11.3 E35-E30 N36-N39
      3 1036 1344 (N20W08 L=122)   367    163        1 ( 958/1317-1125/1362) -21.8 W02-W12 N16-N22
      4 1590 1285 (N05W40 L=155)   669    165        4 (1540/1205-1626/1358) -83.7 W38-W42 N00-N09
      5 1825 1208 (S05W58 L=173)  2574    527        4 (1754/1031-1906/1356)  25.8 W47-W78 S13-N01
      6  182 1029 (N19E59 L= 56)  2277    402        6 ( 115/ 929- 272/1147) -49.1 E65-E51 N11-N26
      7  982  503 (S27E15 L=100)   708    199        1 ( 945/ 428-1013/ 611)  54.0 E18-E10 S32-S20
      8  487  418 (S22E52 L= 63)  1943    615        5 ( 458/ 278- 523/ 600)  39.8 E64-E40 S31-S11
      9  871  233 (S44E36 L= 79)   397     98        0 ( 857/ 208- 888/ 267)  51.3 E38-E33 S46-S41
  

    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