Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  545 1469 (N15E38 L=136)   704    155        1 ( 487/1435- 592/1519)  45.0 E41-E34 N12-N18
      2 1452 1421 (N25W22 L=196)   686     82        7 (1412/1389-1492/1466) -38.7 W19-W25 N22-N27
      3  632 1367 (N10E29 L=145)   236    101        1 ( 618/1331- 654/1436) -71.6 E30-E26 N07-N13
      4 1530 1266 (N16W28 L=202)   400     93        2 (1486/1238-1570/1297) -18.4 W25-W31 N14-N17
      5 1329 1128 (N04W17 L=191)   197     79        0 (1300/1091-1356/1150) -36.9 W14-W18 N02-N05
      6  133  952 (S18E70 L=104)  1022     85        2 ( 109/ 910- 157/ 994)  45.0 E75-E65 S20-S15
      7 1709  809 (S07W49 L=223)   472     40        2 (1680/ 784-1739/ 831)  21.8 W46-W51 S08-S06
      8 1044  737 (S23W05 L=179)   453    169        3 ( 976/ 710-1120/ 765)  18.4 W00-W09 S24-S19
      9 1595  656 (S18W45 L=218)   730    255        0 (1485/ 610-1707/ 688)  31.0 W37-W52 S23-S14
  

    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