Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  836 1510 (N41E01 L=254)   283     60        0 ( 802/1496- 873/1524)   0.0 E03-W02 N39-N41
      2  924 1440 (N34W03 L=259)   506     98        0 ( 884/1390- 964/1482) -74.1 W01-W05 N30-N37
      3  376 1124 (N25E42 L=214)  3039    726        6 ( 183/1005- 598/1257) -31.7 E60-E26 N14-N35
      4 1703  875 (S18W41 L=296)  1989    345        3 (1613/ 748-1755/1048)  43.2 W31-W48 S24-S08
      5  675  975 (N11E22 L=233)   249     58        2 ( 642/ 950- 702/ 996)   0.0 E24-E19 N09-N12
      6  804  934 (N06E15 L=241)   396     63        1 ( 767/ 907- 839/ 962)   0.0 E17-E12 N04-N08
      7  789  733 (S05E20 L=235)   471    162        0 ( 710/ 700- 872/ 763)  -5.2 E25-E14 S06-S04
      8 1182  595 (S22W00 L=255)   250     45        0 (1158/ 566-1201/ 620) -90.0 E01-W01 S24-S20
      9 1532  491 (S40W24 L=279)   415     87        0 (1496/ 464-1570/ 513)  -8.1 W20-W27 S41-S37
     10  546  392 (S23E49 L=206)   460     97        1 ( 513/ 352- 589/ 432) -90.0 E51-E47 S26-S19
  

    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