Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1552 1559 (N41W21 L= 85)   787    100        1 (1510/1516-1593/1608) -33.7 W16-W24 N38-N44
      2  252 1463 (N02E69 L=356)   805     76        1 ( 223/1437- 294/1491)  -8.1 E76-E62 N00-N03
      3 1043 1412 (N15E09 L= 55)   340    104        1 (1011/1359-1063/1457) -39.8 E12-E06 N12-N17
      4 1255 1214 (N09W08 L= 72)   190     96        0 (1244/1166-1272/1256) -53.1 W06-W09 N07-N11
      5 1659 1207 (N21W34 L= 99)   339     43        1 (1634/1183-1677/1237) -56.3 W32-W35 N19-N23
      6  231 1131 (S19E58 L=  7)  1146    210        0 ( 168/1097- 307/1175)  10.6 E66-E50 S20-S15
      7 1869 1099 (N23W57 L=121)   387     63        1 (1844/1062-1891/1126)  56.3 W54-W59 N19-N25
      8 1376 1040 (N03W19 L= 84)   261     63        0 (1351/1007-1407/1070)  63.4 W17-W20 N00-N05
      9  192  874 (S35E62 L=  3)   862    108        1 ( 138/ 849- 244/ 903)   0.0 E71-E54 S35-S32
  

    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