Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1658 1604 (N32W53 L=190)  1984    315        4 (1560/1548-1743/1643) -10.0 W44-W62 N27-N34
      2  800 1595 (N35E13 L=124)   710    165        2 ( 721/1566- 875/1620)   0.0 E18-E07 N33-N37
      3  394 1231 (N14E41 L= 96)   362     92        0 ( 371/1184- 413/1272) -59.0 E42-E38 N11-N16
      4  159 1199 (N14E64 L= 73)   527    143        0 ( 124/1184- 192/1213) -12.5 E68-E59 N12-N15
      5 1028 1132 (N04W01 L=138)   330     82        1 ( 998/1073-1050/1174) -80.5 W00-W01 N00-N06
      6  334 1069 (N05E45 L= 91)   338     78        0 ( 314/1026- 357/1111) -59.0 E46-E43 N02-N07
      7 1369  964 (S07W21 L=158)   201     85        0 (1336/ 944-1410/ 999) -31.0 W18-W23 S08-S05
      8  388  841 (S09E43 L= 94)   292     97        0 ( 348/ 824- 425/ 860)   0.0 E46-E39 S10-S08
      9  290  641 (S21E56 L= 81)  2105    355        7 ( 248/ 550- 342/ 759) -85.6 E60-E52 S26-S13
     10 1761  687 (S26W52 L=189)   772     94        3 (1728/ 653-1796/ 721)  18.4 W49-W55 S27-S23
     11  365  606 (S23E50 L= 87)   347     73        0 ( 339/ 562- 382/ 639) -78.7 E51-E48 S26-S21
     12  904  471 (S36E12 L=124)   229     51        0 ( 872/ 459- 935/ 485)   0.0 E14-E10 S37-S35
  

    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