Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1611 1541 (N21W55 L=133)  2159    415        5 (1454/1502-1750/1604) -36.9 W45-W66 N14-N29
      2  318 1427 (N46E49 L= 29)   815    180        2 ( 281/1389- 355/1469) -28.6 E54-E43 N43-N49
      3 1292 1290 (N16W22 L=100)   499    165        2 (1214/1258-1336/1345)   6.3 W17-W26 N14-N18
      4  327 1133 (N27E46 L= 32)   679    106        0 ( 288/1100- 367/1174) -39.8 E49-E42 N24-N29
      5  911  618 (S14E16 L= 62)   954    218        2 ( 860/ 522- 959/ 703)  48.0 E21-E11 S19-S09
      6  325  554 (S09E62 L= 16)  2411    301        6 ( 254/ 491- 411/ 625) -55.0 E68-E53 S13-S03
      7 1726  588 (S39W42 L=120)   598    155        2 (1679/ 563-1756/ 606) -33.7 W36-W45 S42-S36
      8 1719  536 (S43W42 L=120)   926    176        2 (1677/ 509-1757/ 569) -35.0 W37-W47 S47-S39
      9  520  362 (S24E53 L= 25)   395     92        0 ( 504/ 327- 537/ 395)  68.2 E54-E51 S26-S21
     10  654  314 (S30E45 L= 33)   603    175        1 ( 618/ 280- 674/ 367)  41.2 E50-E41 S32-S25
  

    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