Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  601 1563 (N38E35 L=102)  2436    378        4 ( 444/1463- 784/1627) -13.7 E54-E17 N30-N42
      2  355 1479 (N31E56 L= 81)  1358    213        3 ( 308/1414- 413/1560) -20.0 E70-E47 N26-N35
      3 1152 1432 (N29W09 L=146)  1062    262        2 (1042/1400-1275/1467)  -3.6 W01-W17 N26-N31
      4 1070 1100 (N08W03 L=140)   227     55        0 (1042/1078-1098/1121) -45.0 W01-W04 N06-N09
      5  476  940 (S03E36 L=102)   560    162        0 ( 409/ 902- 527/ 999)  32.0 E40-E31 S05-N01
      6  985  787 (S11E02 L=135)   206     33        0 ( 970/ 758-1003/ 809)  71.6 E03-E01 S12-S09
      7 1068  746 (S14W03 L=140)   299     56        0 (1051/ 714-1092/ 773)  71.6 W01-W04 S15-S11
      8 1613  566 (S27W45 L=182)  2130    363        1 (1558/ 482-1697/ 685) -51.8 W38-W53 S32-S18
      9  901  614 (S23E08 L=129)   471    127        0 ( 840/ 598- 969/ 632)   6.3 E12-E03 S23-S21
     10  644  396 (S39E31 L=106)  1333    227        2 ( 572/ 345- 722/ 488)  32.3 E40-E21 S43-S31
     11  933  215 (S56E09 L=128)   443     37        0 ( 906/ 201- 964/ 231)  -9.5 E12-E06 S57-S54
  

    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