Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1326 1771 (N55W36 L=133)   526     52        1 (1298/1750-1356/1791)  16.7 W31-W41 N53-N57
      2  856 1440 (N30E11 L= 86)  1236    294        2 ( 757/1359- 973/1520) -34.5 E19-E03 N24-N35
      3 1606 1456 (N29W46 L=143)  1773    302        1 (1462/1429-1736/1479)   0.0 W32-W61 N27-N30
      4  211 1340 (N22E68 L= 30)  2224    174        7 ( 154/1291- 260/1409)  82.9 E80-E59 N18-N26
      5  202 1047 (N04E61 L= 36)   483     85        0 ( 165/1023- 244/1060)  11.3 E66-E56 N02-N05
      6  549  932 (S02E30 L= 67)   227     43        0 ( 536/ 900- 561/ 959) -71.6 E31-E29 S03-N00
      7 1679  855 (S08W44 L=142)   332     63        0 (1662/ 816-1696/ 883) -53.1 W42-W45 S10-S06
      8  808  714 (S15E14 L= 83)   434    118        1 ( 782/ 661- 842/ 768)  60.3 E16-E11 S18-S11
      9 1637  696 (S18W43 L=140)   290     49        1 (1619/ 673-1665/ 717)  14.0 W41-W45 S19-S16
     10 1192  392 (S39W12 L=110)   583    126        2 (1132/ 353-1249/ 435)  29.1 W07-W16 S41-S35
     11 1265  198 (S59W28 L=125)   533     46        0 (1237/ 181-1289/ 214) -18.4 W23-W32 S61-S56
  

    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