Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1107 1904 (N59E39 L=194)   577     52        1 (1085/1882-1125/1928) -38.7 E45-E34 N55-N63
      2 1176 1858 (N57E25 L=208)   873    101        0 (1142/1817-1207/1886) -11.3 E30-E19 N52-N60
      3  838 1674 (N29E33 L=200)  1881    484        3 ( 715/1538- 970/1821) -18.4 E54-E18 N23-N37
      4 1410 1756 (N58W04 L=237)   418     53        0 (1387/1733-1439/1772)  50.2 W02-W07 N54-N60
      5 1292 1509 (N32W03 L=236)   789    127        4 (1244/1456-1339/1564) -29.7 E01-W06 N29-N34
      6  964 1415 (N16E14 L=219)   302    117        1 ( 916/1373-1014/1446)   0.0 E17-E10 N15-N17
      7  226 1283 (S08E64 L=169)  1144    115        1 ( 189/1231- 264/1329)  45.0 E67-E60 S11-S05
      8 1264 1244 (N14W08 L=241)   551    156        2 (1240/1165-1297/1314) -78.7 W06-W09 N09-N19
      9  577  990 (S18E26 L=207)   404    106        1 ( 519/ 971- 628/1001)  14.0 E30-E22 S18-S15
     10 1007  780 (S18W06 L=239)   362    124        0 ( 954/ 760-1067/ 807)   7.1 W01-W09 S19-S17
     11 1441  644 (S14W36 L=269)   320    118        0 (1415/ 591-1458/ 690) -39.8 W32-W38 S16-S11
     12 1713  649 (S05W57 L=290)   359     30        0 (1696/ 628-1733/ 670) -11.3 W54-W59 S06-S04
  

    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