Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1166 1677 (N38E02 L=244)   934    221        0 (1064/1625-1263/1721) -11.3 E10-W05 N35-N39
      2 1007 1537 (N25E10 L=236)   447    185        0 ( 934/1483-1079/1578) -20.0 E15-E04 N22-N26
      3 1706 1448 (N33W43 L=290)   640    146        0 (1665/1389-1745/1513) -53.1 W40-W47 N29-N37
      4  751 1075 (S07E17 L=229)   191     44        1 ( 733/1047- 772/1098)  90.0 E17-E15 S09-S05
      5  342  904 (S22E45 L=201)   802    197        0 ( 296/ 819- 374/ 988)  68.2 E50-E42 S27-S17
      6 1011  931 (S12W01 L=247)   203    144        0 ( 951/ 887-1069/ 966)  54.5 E02-W03 S15-S08
      7 1335  805 (S14W22 L=269)   525    165        0 (1264/ 768-1413/ 839)  41.2 W18-W26 S17-S10
      8  923  665 (S29E00 L=246)  2252    661        4 ( 622/ 587-1207/ 724)  20.7 E21-W16 S36-S20
      9 1710  651 (S15W55 L=301)   454     48        1 (1686/ 630-1732/ 677)  71.6 W52-W57 S16-S13
     10 1644  558 (S22W54 L=300)   754    136        1 (1574/ 539-1707/ 583)  24.4 W47-W58 S23-S18
     11 1046  201 (S61W30 L=277)  1093    114        1 ( 992/ 154-1080/ 253) -17.7 W19-W41 S65-S57
  

    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