Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1477 1734 (N58W39 L=172)  1408    213        3 (1388/1706-1578/1754)   3.5 W25-W60 N53-N60
      2  901 1642 (N40E19 L=114)   409     95        1 ( 857/1608- 929/1662) -23.2 E23-E16 N37-N40
      3 1411 1626 (N46W24 L=158)   322     36        0 (1385/1611-1436/1639) -14.0 W21-W27 N45-N47
      4  645 1450 (N23E32 L=102)   731    164        2 ( 602/1396- 681/1526) -83.7 E34-E28 N18-N27
      5 1196 1458 (N30W06 L=140)   272     87        0 (1182/1423-1213/1505) -59.0 W04-W07 N28-N33
      6  475 1427 (N19E44 L= 89)   851    162        5 ( 427/1395- 532/1468) -10.3 E49-E38 N16-N21
      7  543 1365 (N16E37 L= 97)   736    217        2 ( 497/1269- 583/1452)  85.2 E38-E35 N09-N21
      8 1723 1421 (N35W53 L=186)   826    117        0 (1663/1393-1774/1449)  -9.5 W46-W59 N33-N35
      9  439  818 (S19E37 L= 96)   954    146        3 ( 395/ 753- 498/ 863)  37.9 E41-E32 S23-S15
     10 1697  651 (S14W54 L=188)   429     30        0 (1673/ 629-1722/ 672)  26.6 W52-W56 S15-S12
     11 1360  392 (S36W37 L=170)   362     93        0 (1319/ 363-1400/ 424)  51.3 W34-W38 S38-S33
  

    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