Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  461 1445 (N46E34 L= 54)   576    181        2 ( 378/1419- 543/1468)   0.0 E43-E25 N44-N46
      2  999 1040 (N09E01 L= 87)   266     84        1 ( 975/1004-1018/1085)  45.0 E03-W00 N06-N11
      3  278  921 (N15E54 L= 35)  1398    227        4 ( 232/ 834- 329/1028) -65.2 E57-E50 N09-N22
      4  191  891 (N15E65 L= 23)  1967    269        4 ( 131/ 820- 240/ 994) -79.7 E75-E59 N09-N20
      5 1498  831 (S16W24 L=112)   808    199        1 (1446/ 747-1555/ 909)  35.0 W19-W29 S19-S12
      6 1106  633 (S17E04 L= 84)   998    432        1 (1048/ 484-1190/ 791)  50.5 E11-W04 S26-S09
      7 1303  606 (S24W07 L= 95)   830    222        1 (1222/ 548-1398/ 672)  12.1 W00-W14 S26-S22
      8 1704  645 (S35W39 L=128)   376     19        0 (1686/ 625-1721/ 670)  71.6 W37-W41 S36-S32
      9  275  614 (S03E65 L= 24)   555     76        1 ( 257/ 582- 301/ 651) -78.7 E67-E62 S05-N00
     10 1437  513 (S35W14 L=103)   509    120        1 (1391/ 477-1481/ 555)  20.6 W10-W18 S36-S32
     11 1133  456 (S30E08 L= 81)   448     54        0 (1096/ 437-1163/ 476)   0.0 E10-E05 S31-S28
     12  735  273 (S34E45 L= 44)  1971    288        1 ( 681/ 163- 785/ 412)  43.7 E62-E38 S44-S22
  

    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