Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  225 1386 (N25E59 L= 67)   505    143        0 ( 192/1369- 255/1404) -14.0 E63-E55 N23-N26
      2  518 1285 (N16E31 L= 94)   362     97        0 ( 471/1258- 554/1303) -26.6 E34-E28 N14-N17
      3  170 1285 (N19E63 L= 62)   502    114        0 ( 143/1269- 200/1300)   0.0 E67-E59 N17-N20
      4  229 1228 (N15E56 L= 70)   437     83        0 ( 203/1203- 254/1254) -31.0 E58-E53 N13-N16
      5 1208 1160 (N05W12 L=137)   256     73        0 (1171/1123-1237/1190) -53.1 W09-W13 N02-N06
      6 1691  955 (S09W44 L=169)   998    238        0 (1655/ 831-1716/1074) -78.7 W41-W46 S16-S01
      7 1514 1017 (S05W30 L=156)   249     56        0 (1495/ 983-1532/1051) -63.4 W29-W31 S06-S02
      8  124  945 (S01E70 L= 55)   726    211        1 (  84/ 922- 163/ 968)  -4.1 E78-E64 S02-N00
      9  497  655 (S22E38 L= 88)   488     94        0 ( 475/ 605- 519/ 700)  80.5 E38-E36 S24-S18
     10 1843  667 (S27W62 L=187)   557     27        0 (1825/ 645-1862/ 686)  56.3 W59-W63 S28-S25
     11  382  606 (S24E49 L= 77)   305     93        0 ( 359/ 580- 421/ 647) -53.1 E49-E45 S25-S21
     12 1089  491 (S36W02 L=127)   656    140        1 (1021/ 467-1168/ 508)   5.7 E03-W07 S37-S34
  

    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