Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1002 1722 (N46E17 L=145)   341     55        0 ( 970/1709-1035/1740)   0.0 E19-E14 N43-N47
      2  355 1464 (N16E56 L=105)  2543    467        8 ( 290/1309- 408/1615) -76.0 E61-E52 N06-N26
      3  558 1301 (N09E34 L=127)   629    181        1 ( 512/1260- 633/1343) -11.3 E38-E28 N06-N11
      4 1090 1253 (N15W00 L=162)   344     83        0 (1045/1218-1127/1282) -21.8 E02-W03 N12-N16
      5  667 1184 (N04E24 L=137)   218     69        0 ( 650/1148- 683/1215) -45.0 E25-E22 N02-N05
      6  585  955 (S11E26 L=135)  1000    245        2 ( 550/ 851- 633/1081) -65.2 E29-E22 S16-S03
      7 1274  955 (S01W16 L=178)   464    126        0 (1211/ 923-1335/ 995)  35.5 W12-W19 S03-N02
      8 1794  973 (N09W55 L=216)   325     54        1 (1773/ 954-1813/ 993)  26.6 W52-W56 N07-N09
      9 1406  591 (S21W33 L=194)   693    141        1 (1322/ 524-1535/ 684)  47.5 W28-W39 S25-S13
     10  233  634 (S38E58 L=103)   717    160        0 ( 208/ 608- 253/ 668)  35.5 E62-E55 S40-S35
     11 1288  416 (S34W31 L=193)   344     81        0 (1259/ 395-1325/ 442) -15.9 W27-W34 S35-S32
     12 1564  371 (S33W56 L=217)  1289    220        0 (1498/ 337-1665/ 428)  32.5 W49-W62 S35-S26
  

    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