Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1232 1674 (N39E08 L=339)   320     72        0 (1203/1658-1254/1701)  63.4 E10-E06 N37-N41
      2  351 1631 (N13E69 L=278)   985     93        2 ( 319/1597- 383/1666)  78.7 E72-E66 N10-N15
      3  750 1637 (N22E36 L=312)   285     61        1 ( 725/1606- 774/1655)  18.4 E37-E34 N19-N22
      4  274 1285 (S10E57 L=290)   345     58        1 ( 244/1257- 297/1316)  53.1 E59-E55 S11-S07
      5  387 1174 (S14E43 L=304)   419    107        2 ( 371/1113- 400/1230) -66.8 E45-E41 S17-S10
      6 1317 1138 (N08W13 L=  0)   255     87        0 (1299/1093-1336/1187) -80.5 W12-W14 N05-N11
      7 1415 1098 (N09W20 L=  7)   218     59        2 (1394/1073-1441/1123) -56.3 W18-W21 N07-N10
      8 1854 1076 (N23W54 L= 41)   593     87        1 (1821/1045-1882/1108)  -9.5 W49-W57 N21-N24
      9  387 1029 (S22E40 L=307)   909    255        2 ( 298/1002- 485/1072)  23.2 E48-E33 S24-S17
     10 1114  518 (S31W22 L= 10)  1558    428        4 ( 940/ 440-1291/ 612)   3.9 W07-W36 S33-S28
  

    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