Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  700 1706 (N48E33 L=342)   507    103        1 ( 676/1670- 721/1742) -63.4 E35-E29 N45-N51
      2  655 1616 (N41E33 L=343)   357     48        1 ( 634/1597- 676/1637)  45.0 E34-E30 N39-N42
      3 1013 1488 (N33E02 L= 13)   642    126        0 ( 952/1451-1080/1522) -14.0 E06-W02 N30-N35
      4 1295 1414 (N28W18 L= 34)   284     47        1 (1277/1387-1317/1441) -63.4 W16-W19 N26-N30
      5 1547 1322 (N22W36 L= 51)  1822    402        4 (1377/1237-1695/1416) -25.6 W23-W47 N16-N28
      6  586 1347 (N22E31 L=345)   412     77        0 ( 567/1311- 607/1392) -78.7 E32-E29 N19-N24
      7 1123 1261 (N18W06 L= 21)   985    283        5 (1065/1130-1175/1371) -68.2 W01-W09 N09-N24
      8 1284 1255 (N18W16 L= 32)   532    162        0 (1221/1189-1344/1339) -52.1 W12-W20 N13-N23
      9  722 1289 (N18E21 L=355)   245     56        1 ( 699/1262- 758/1307)  45.0 E22-E18 N16-N19
     10  924 1066 (N05E06 L=  9)  1331    356        1 ( 838/ 973-1072/1195) -37.6 E12-W02 N00-N13
     11  771  921 (S04E15 L=  0)   354    153        1 ( 697/ 908- 844/ 937)   5.7 E20-E10 S04-S02
     12 1103  246 (S53W11 L= 27)   373     90        0 (1054/ 232-1146/ 256)  -5.7 W06-W16 S54-S52
  

    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