Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  360 1556 (N56E56 L=101)   666     65        0 ( 329/1531- 387/1584)  36.9 E62-E51 N54-N57
      2 1783 1374 (N07W64 L=220)   744     96        1 (1753/1330-1814/1410)  11.3 W58-W68 N05-N09
      3  516 1373 (N39E29 L=128)   295     45        0 ( 493/1352- 542/1394)  11.3 E31-E26 N37-N40
      4  985 1360 (N28W06 L=162)   370     87        1 ( 942/1334-1035/1377)   0.0 W02-W09 N26-N28
      5 1788 1035 (S12W51 L=208)   671     88        1 (1758/ 991-1816/1071) -90.0 W48-W55 S14-S09
      6 1711  807 (S23W41 L=198)   307     68        0 (1674/ 777-1737/ 833)   9.5 W37-W43 S24-S21
      7  340  517 (S11E63 L= 93)   665     71        2 ( 308/ 493- 363/ 542) -21.8 E68-E60 S12-S09
      8  513  417 (S20E52 L=105)  1752    258        1 ( 457/ 315- 578/ 512)  24.0 E69-E42 S27-S14
      9 1043  446 (S29E13 L=144)   276    126        0 ( 987/ 431-1113/ 461)  -7.1 E16-E08 S29-S27
     10  606  362 (S26E48 L=109)   553    116        1 ( 567/ 341- 667/ 380) -23.2 E51-E43 S27-S24
     11  731  318 (S31E41 L=116)   378     65        0 ( 700/ 297- 759/ 342)   8.1 E44-E37 S32-S29
  

    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