Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  355 1559 (N56E57 L=100)   560     73        0 ( 321/1538- 382/1584)  12.5 E64-E52 N54-N58
      2 1776 1367 (N07W62 L=220)   708     90        0 (1745/1324-1806/1398)  10.3 W57-W68 N05-N08
      3  509 1372 (N39E30 L=128)   251     41        0 ( 488/1353- 530/1393)  14.0 E31-E27 N37-N40
      4  977 1355 (N28W05 L=162)   316     93        1 ( 923/1337-1019/1373)  -8.1 W01-W08 N25-N28
      5 1776 1037 (S11W50 L=208)   562     74        0 (1747/1006-1801/1074)  21.8 W47-W53 S12-S09
      6  339  518 (S11E64 L= 94)   513     45        0 ( 321/ 491- 357/ 545)  76.0 E66-E61 S13-S09
      7  505  417 (S20E53 L=105)  1652    262        1 ( 453/ 313- 567/ 513)  23.2 E71-E43 S27-S13
      8 1031  443 (S29E14 L=144)   253    126        0 ( 975/ 429-1097/ 457)  -7.1 E17-E09 S29-S27
      9  594  363 (S25E49 L=109)   581    111        2 ( 560/ 341- 657/ 382) -33.7 E53-E44 S27-S23
     10  724  318 (S31E41 L=116)   334     65        0 ( 693/ 299- 749/ 343)  15.9 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