Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1541 1693 (N32W58 L= 26)   652    124        1 (1510/1662-1579/1712) -26.6 W54-W61 N30-N33
      2  985 1571 (N41W12 L=340)   384     63        2 ( 953/1540-1009/1600)   9.5 W08-W14 N38-N42
      3 1734 1416 (N12W60 L= 27)   401    110        1 (1717/1384-1750/1444)  63.4 W57-W62 N09-N13
      4  991 1404 (N30W07 L=335)   204     20        0 ( 968/1390-1011/1421)  18.4 W05-W08 N29-N31
      5 1424 1399 (N19W34 L=  2)   339     83        1 (1383/1381-1463/1417) -18.4 W31-W37 N17-N20
      6 1254 1240 (N14W18 L=346)   276     74        0 (1210/1207-1283/1266) -53.1 W16-W19 N12-N16
      7  252  850 (N10E58 L=270)   439     89        0 ( 226/ 818- 280/ 899) -59.0 E59-E55 N07-N12
      8 1404  713 (S20W16 L=344)   385    117        1 (1370/ 657-1445/ 765)  45.0 W13-W19 S22-S17
      9 1212  426 (S35E02 L=326)   317     45        0 (1189/ 401-1239/ 449)  45.0 E03-W00 S36-S33
     10 1510  336 (S53W23 L=351)   411     71        1 (1477/ 322-1543/ 354)  -8.1 W19-W27 S55-S51
  

    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