Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  504 1495 (N31E43 L= 37)   662    124        1 ( 446/1473- 565/1513)  -4.8 E48-E36 N29-N32
      2  917 1362 (N24E08 L= 71)   467    163        0 ( 850/1351-1009/1375)   0.0 E12-E01 N23-N24
      3 1481 1202 (N15W30 L=109)   297    129        0 (1422/1184-1522/1228) -23.2 W25-W32 N13-N16
      4  160  869 (S11E69 L= 10)  2146    260        7 ( 120/ 781- 221/ 955) -39.8 E76-E62 S15-S05
      5 1409  840 (S08W25 L=104)   271     53        1 (1383/ 812-1434/ 863)  45.0 W23-W26 S09-S06
      6 1652  734 (S14W45 L=124)  1094    187        3 (1593/ 671-1717/ 817)  41.6 W40-W49 S18-S09
      7  372  701 (S20E47 L= 33)   998    185        2 ( 340/ 616- 423/ 786) -79.7 E49-E43 S25-S14
      8 1063  496 (S31W04 L= 84)   556    138        1 (1041/ 429-1088/ 558) -78.7 W02-W06 S36-S26
      9 1111  375 (S41W09 L= 89)   336     49        0 (1088/ 351-1136/ 405)  53.1 W07-W10 S42-S38
  

    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