Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1413 1453 (N18W36 L=327)   280     77        0 (1379/1429-1438/1486)   9.5 W33-W39 N16-N19
      2 1109 1369 (N22W15 L=305)   230     54        0 (1085/1342-1142/1390)   0.0 W12-W17 N21-N23
      3 1191 1338 (N18W18 L=309)   304     67        0 (1154/1314-1232/1366) -53.1 W17-W20 N15-N20
      4 1297 1131 (N04W18 L=308)   405    107        2 (1257/1088-1355/1170)  15.9 W15-W22 N02-N05
      5 1169 1146 (N08W11 L=302)   494    121        0 (1094/1128-1234/1163) -20.6 W06-W14 N05-N10
      6  828 1073 (N13E10 L=281)   537    179        1 ( 732/1019- 911/1112)   0.0 E16-E04 N12-N14
      7  188  877 (N16E62 L=228)  2617    490        2 (  88/ 829- 299/ 924) -14.0 E75-E51 N12-N19
      8 1339  650 (S24W08 L=299)  1581    368        6 (1286/ 496-1400/ 792)  45.0 W00-W15 S33-S17
      9  580  569 (S09E40 L=251)   257     57        0 ( 558/ 542- 597/ 595) -71.6 E40-E38 S10-S07
     10  735  508 (S17E32 L=258)   276     40        0 ( 707/ 492- 762/ 523) -33.7 E33-E30 S17-S15
     11  587  350 (S23E50 L=240)   427     36        1 ( 558/ 332- 610/ 365) -56.3 E52-E47 S24-S21
  

    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