Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  466 1579 (N15E50 L=227)   428     50        1 ( 443/1554- 489/1604)  90.0 E52-E48 N13-N17
      2 1190 1423 (N20W01 L=278)  1027    201        2 (1116/1338-1253/1495) -31.0 E04-W06 N15-N23
      3 1197 1193 (N06W06 L=284)   358     96        1 (1161/1155-1231/1235) -31.0 W03-W08 N04-N08
      4 1742  984 (N08W46 L=324)   781    216        0 (1708/ 899-1805/1022)  90.0 W43-W52 N02-N10
      5  240  902 (S28E52 L=225)   358     76        0 ( 211/ 878- 269/ 928)  31.0 E55-E49 S29-S26
      6  294  821 (S33E47 L=230)   388     96        0 ( 258/ 807- 332/ 836)   8.1 E50-E43 S33-S31
      7 1719  550 (S15W61 L=338)   457     47        0 (1692/ 517-1739/ 580)  76.0 W59-W62 S17-S13
      8 1466  478 (S27W44 L=321)  1104    229        1 (1357/ 446-1576/ 515)  34.7 W37-W50 S31-S22
      9  863  208 (S66W20 L=297)   570     77        1 ( 840/ 181- 887/ 237) -21.8 W15-W25 S68-S62
  

    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