Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1421 1684 (N30W51 L= 79)   736     90        1 (1382/1654-1467/1714) -68.2 W48-W54 N27-N32
      2  536 1555 (N51E21 L=  7)   352     49        0 ( 514/1531- 566/1572)   0.0 E23-E17 N49-N53
      3  201 1307 (N42E63 L=326)   728     77        0 ( 156/1285- 241/1323)  -6.3 E73-E55 N40-N43
      4 1201 1223 (N11W15 L= 44)   301    117        0 (1173/1173-1227/1272)  80.5 W14-W16 N08-N14
      5  514 1157 (N25E28 L=  0)   426    127        2 ( 459/1119- 572/1175) -36.9 E32-E24 N21-N27
      6  170 1112 (N30E64 L=325)  1221    127        1 ( 120/1049- 220/1165) -29.4 E74-E56 N25-N34
      7  325  997 (N20E46 L=342)   476    103        0 ( 265/ 983- 373/1014) -10.3 E52-E41 N19-N21
      8 1724  959 (S18W42 L= 70)   289     35        0 (1707/ 935-1745/ 979)  14.0 W40-W44 S19-S17
      9 1014  756 (S09E08 L= 21)   319    122        0 ( 951/ 743-1072/ 769) -15.9 E11-E04 S10-S07
     10 1287  430 (S38E01 L= 27)  2324    523        3 (1063/ 315-1459/ 579)   8.1 E21-W14 S41-S32
     11  446  496 (S11E54 L=335)   393     50        1 ( 414/ 477- 470/ 515) -21.8 E55-E50 S12-S09
     12  878  353 (S31E31 L=358)  1707    285        2 ( 796/ 260- 950/ 463)  24.4 E43-E21 S35-S24
  

    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