Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  851 1840 (N55E19 L=162)   615     80        1 ( 806/1812- 892/1857)  12.5 E23-E14 N51-N57
      2 1632 1509 (N28W46 L=227)   356     65        0 (1595/1485-1661/1528) -31.0 W42-W48 N26-N29
      3  297 1453 (N25E57 L=124)  1081    110        2 ( 267/1400- 326/1510) -48.8 E61-E51 N20-N28
      4  462 1415 (N21E39 L=142)   464    137        0 ( 424/1366- 493/1468) -41.2 E43-E35 N18-N25
      5  943 1322 (N15E05 L=176)   289     83        0 ( 902/1305- 993/1340) -18.4 E07-E01 N13-N15
      6  388 1241 (N10E42 L=139)   360     56        2 ( 373/1200- 407/1276)  90.0 E43-E40 N07-N12
      7  463 1121 (N03E36 L=145)  1569    286        6 ( 391/ 980- 543/1229)  54.5 E41-E30 S05-N09
      8 1119  834 (S15W06 L=187)  1054    349        1 (1036/ 697-1173/ 978)  64.8 W00-W09 S23-S06
      9 1629  786 (S17W41 L=222)   395    139        0 (1575/ 765-1678/ 813)   0.0 W36-W45 S17-S15
     10 1092  473 (S38W05 L=186)  1237    406        2 ( 917/ 446-1299/ 505)   8.1 E08-W20 S40-S35
     11 1004  305 (S52E02 L=179)   368     36        0 ( 975/ 287-1036/ 322)  11.3 E04-W01 S53-S50
  

    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