Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  379 1532 (N30E51 L=268)  3069    346        2 ( 316/1385- 428/1668) -30.1 E81-E42 N19-N41
      2 1465 1289 (N13W28 L=347)   250     64        0 (1440/1264-1495/1316) -45.0 W26-W29 N11-N14
      3 1435 1135 (N04W25 L=344)   501    111        2 (1401/1082-1471/1184)  71.6 W23-W27 N00-N06
      4  778  901 (S11E15 L=304)   258     91        0 ( 730/ 882- 826/ 919) -21.8 E18-E12 S11-S09
      5  413  834 (S14E41 L=278)   578    159        0 ( 384/ 765- 470/ 881)  36.9 E44-E36 S18-S11
      6  792  820 (S16E14 L=305)   291     55        1 ( 762/ 793- 820/ 846) -36.9 E16-E12 S17-S13
      7 1716  649 (S25W53 L= 12)   308     40        0 (1703/ 623-1735/ 676) -36.9 W50-W55 S26-S23
      8 1107  442 (S41W07 L=326)   392    113        1 (1048/ 421-1166/ 455)   6.3 W02-W11 S42-S39
      9  935  225 (S60E10 L=309)  1397    200        1 ( 844/ 191-1032/ 267)  12.3 E22-W01 S63-S55
  

    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