Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1063 1391 (N14E07 L=238)  1289    300        1 (1009/1264-1118/1537) -47.1 E13-W00 N07-N22
      2  450 1409 (N02E46 L=200)  1136    283        2 ( 326/1382- 553/1427)  15.5 E55-E37 S01-N04
      3 1678 1328 (N29W34 L=280)   343    103        0 (1630/1306-1708/1341)  26.6 W30-W36 N26-N30
      4 1630 1277 (N24W31 L=276)   275     32        0 (1610/1254-1648/1301) -56.3 W29-W32 N22-N25
      5  573 1256 (S04E32 L=213)   213     18        1 ( 558/1233- 594/1277) -45.0 E33-E30 S05-S03
      6 1204 1183 (N06W06 L=252)   673    140        2 (1125/1126-1269/1232)  -5.7 W00-W10 N03-N07
      7  471 1161 (S12E37 L=209)   479    108        2 ( 418/1130- 522/1189)  29.7 E40-E33 S14-S09
      8 1638  990 (N07W37 L=283)  1459    264        1 (1599/ 871-1678/1125) -42.5 W31-W43 N01-N13
      9  268 1084 (S20E52 L=194)   388    136        0 ( 236/1061- 312/1121) -23.2 E55-E48 S20-S17
     10  840  807 (S24E05 L=240)  1144    399        0 ( 649/ 770-1031/ 852)  26.6 E17-W05 S29-S18
     11  339  756 (S39E42 L=203)  1587    396        3 ( 224/ 705- 458/ 788)  21.4 E54-E31 S43-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