Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  928 1619 (N37E17 L=354)   425    147        0 ( 885/1568- 985/1671) -28.6 E22-E11 N34-N40
      2 1264 1549 (N38W10 L= 21)   282     45        1 (1243/1533-1296/1569) -11.3 W07-W12 N36-N38
      3  502 1378 (N15E40 L=331)   664    202        1 ( 470/1283- 528/1464) -54.0 E45-E36 N09-N20
      4 1732 1446 (N37W54 L= 66)   476     41        0 (1697/1431-1759/1463)   6.3 W49-W58 N36-N38
      5  232 1366 (N11E67 L=304)   932    104        1 ( 204/1330- 251/1415)  78.7 E70-E63 N08-N13
      6 1652 1348 (N29W42 L= 53)   431     94        1 (1615/1311-1689/1384) -23.2 W38-W45 N26-N31
      7  747 1231 (N09E20 L=351)   874    204        1 ( 704/1132- 787/1317)  85.6 E21-E18 N02-N15
      8  970 1178 (N09E05 L=  6)   332    150        0 ( 933/1119-1012/1234) -35.5 E08-E01 N05-N12
      9  431 1115 (S02E40 L=331)   411     77        0 ( 403/1083- 464/1155) -36.9 E42-E37 S03-N00
     10  955  662 (S23W01 L= 12)  2215    662        3 ( 698/ 565-1232/ 767)  32.0 E16-W16 S32-S12
     11  333  596 (S37E52 L=320)   988    208        0 ( 288/ 526- 395/ 690) -55.0 E55-E46 S41-S30
     12  477  632 (S32E35 L=337)   406     77        0 ( 434/ 620- 525/ 648)   0.0 E38-E30 S32-S30
  

    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