Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  721 1562 (N24E28 L=170)   236     83        2 ( 702/1530- 748/1598) -45.0 E29-E24 N21-N26
      2  328 1523 (N18E60 L=139)  2442    160        7 ( 256/1462- 385/1591)  28.6 E65-E54 N13-N22
      3 1089 1514 (N25E03 L=196)   806    141        2 (1048/1447-1117/1571) -53.1 E06-W00 N20-N28
      4  344 1427 (N12E53 L=145)   735     67        2 ( 309/1392- 374/1454) -33.7 E56-E49 N09-N13
      5 1765 1286 (N22W49 L=247)  2359    381        2 (1606/1154-1862/1403) -21.8 W34-W59 N15-N26
      6 1162 1358 (N16W04 L=202)   424    104        1 (1109/1324-1213/1385)   0.0 W00-W07 N14-N17
      7  527 1233 (N01E34 L=165)   389    120        0 ( 487/1208- 575/1262)  18.4 E36-E30 N00-N03
      8  930 1199 (N03E08 L=190)   387    101        1 ( 893/1156- 966/1240) -38.7 E10-E05 N01-N05
      9  390 1122 (S07E42 L=156)   272     89        0 ( 365/1088- 429/1146)  38.7 E44-E39 S08-S04
     10 1367  908 (S08W22 L=220)   244     67        0 (1346/ 882-1393/ 943) -26.6 W20-W24 S09-S06
     11  671  803 (S23E20 L=178)   469    125        0 ( 613/ 785- 733/ 829)   0.0 E23-E15 S24-S21
     12 1307  748 (S18W21 L=220)   705    289        1 (1171/ 699-1426/ 789)  28.1 W13-W28 S22-S14
     13 1692  417 (S30W71 L=269)  2647    235        2 (1598/ 362-1789/ 493)  38.7 W63-W79 S35-S23
  

    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