Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1077 1802 (N44E20 L=353)   515     61        0 (1038/1783-1114/1820)   9.5 E22-E16 N42-N46
      2 1310 1681 (N41W02 L= 15)   343     36        0 (1288/1658-1331/1702) -45.0 W00-W03 N39-N43
      3  585 1568 (N16E41 L=332)   598    167        1 ( 561/1500- 603/1640) -52.1 E44-E37 N11-N20
      4 1479 1577 (N38W17 L= 30)   440     83        0 (1446/1556-1521/1599)   0.0 W14-W20 N35-N40
      5 1185 1335 (N14W02 L= 16)   967    284        1 (1138/1219-1271/1464) -39.8 E03-W09 N08-N21
      6 1887 1160 (N24W57 L= 70)   658     93        1 (1857/1128-1913/1205) -18.4 W53-W60 N22-N25
      7  254 1083 (S18E53 L=320)   946    161        1 ( 222/ 994- 273/1153) -60.9 E57-E50 S23-S14
      8  729 1042 (S12E17 L=356)   904    191        0 ( 647/ 961- 790/1129)  71.6 E21-E14 S18-S06
      9 1447 1098 (N07W23 L= 36)   315    111        0 (1404/1046-1470/1128) -21.8 W19-W24 N04-N08
     10 1351  888 (S07W22 L= 35)   993    241        1 (1267/ 775-1442/1003)  76.0 W19-W24 S15-N01
     11 1027  736 (S23W07 L= 20)   719    163        2 ( 938/ 712-1110/ 766)   9.5 W00-W12 S25-S21
     12 1784  642 (S08W63 L= 76)  1579    222        1 (1707/ 575-1844/ 723)  59.7 W56-W67 S13-S01
  

    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