Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  940 1528 (N37W11 L=144)   959    292        2 ( 805/1492-1077/1545) -27.9 W01-W18 N30-N40
      2 1578 1066 (S07W33 L=166)   244     19        0 (1559/1046-1595/1087)   0.0 W31-W34 S08-S06
      3 1565  903 (S16W28 L=161)   976    366        1 (1431/ 802-1670/1031)  13.4 W17-W38 S18-S12
      4 1197  848 (S08W05 L=137)   971    260        2 (1133/ 726-1248/ 959)  86.4 W02-W06 S16-N00
      5  245  887 (N16E57 L= 76)   614    101        2 ( 210/ 851- 273/ 931) -90.0 E60-E53 N13-N18
      6 1359  842 (S13W14 L=147)   240    107        0 (1331/ 803-1379/ 874)  63.4 W11-W15 S15-S11
      7 1586  800 (S23W28 L=161)   274    121        0 (1556/ 752-1620/ 834)  18.4 W24-W30 S25-S21
      8 1602  705 (S29W27 L=160)   315     45        1 (1575/ 689-1630/ 727) -36.9 W25-W29 S31-S27
      9 1071  641 (S17E08 L=125)   297     64        1 (1041/ 622-1113/ 667) -36.9 E09-E05 S18-S14
     10 1400  580 (S31W10 L=143)   294     37        0 (1373/ 562-1426/ 601)   0.0 W08-W12 S31-S29
     11 1275  461 (S35E01 L=131)  1164    204        1 (1219/ 382-1322/ 534)  52.1 E07-W02 S39-S28
     12  568  448 (S16E46 L= 86)   369     70        0 ( 549/ 414- 584/ 483)  51.3 E48-E43 S18-S13
  

    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