Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  391 1424 (N30E43 L= 16)  3872    495        8 ( 265/1242- 499/1613) -35.9 E77-E31 N16-N45
      2  746 1440 (N27E15 L= 45)   471    196        1 ( 689/1353- 812/1507) -56.3 E18-E10 N20-N32
      3 1556 1486 (N22W42 L=101)   777    191        1 (1457/1456-1635/1507) -13.0 W34-W48 N19-N24
      4  166 1208 (N18E64 L=355)  5328    566        8 (  83/1003- 235/1389)  80.5 E79-E56 N05-N29
      5  223  753 (S10E62 L=358)   556    108        0 ( 205/ 701- 244/ 806)  40.6 E65-E57 S12-S06
      6  439  705 (S15E42 L= 17)   902    274        0 ( 363/ 600- 510/ 803)  34.5 E52-E35 S20-S09
      7 1286  745 (S20W14 L= 73)   345    140        1 (1231/ 714-1355/ 769)  14.0 W10-W18 S21-S19
      8 1092  669 (S23W01 L= 60)   748    249        2 ( 992/ 639-1191/ 732) -32.5 E04-W07 S25-S18
      9  251  561 (S22E69 L=350)  2231    172        4 ( 204/ 466- 296/ 614)  84.3 E75-E64 S27-S17
     10 1473  464 (S41W30 L= 89)   952    101        2 (1428/ 429-1523/ 504)   0.0 W25-W34 S43-S38
     11 1288  415 (S43W14 L= 73)   252     39        1 (1265/ 396-1313/ 441) -36.9 W12-W16 S44-S40
  

    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