Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1004 1880 (N60W47 L=  7)   966    127        2 ( 966/1832-1042/1924)  16.9 W38-W61 N55-N62
      2  838 1766 (N58W17 L=337)  2436    336        4 ( 701/1658- 956/1855)   8.7 E01-W38 N53-N62
      3  380 1533 (N56E43 L=277)   363     47        0 ( 358/1508- 402/1555)  56.3 E46-E39 N53-N57
      4  961 1436 (N31W08 L=328)   360    121        0 ( 914/1400- 986/1484)  23.2 W04-W11 N29-N33
      5  264 1262 (N39E52 L=268)   517     70        3 ( 222/1241- 299/1275) -15.3 E58-E47 N36-N39
      6 1616  878 (S19W32 L=352)   557    124        0 (1599/ 824-1631/ 945)  50.2 W29-W34 S22-S15
      7 1453  750 (S21W18 L=338)   314    134        0 (1420/ 697-1481/ 810)  39.8 W14-W20 S23-S18
      8 1514  501 (S40W17 L=337)  1470    215        3 (1470/ 439-1582/ 588)  59.0 W13-W23 S45-S34
      9 1101  487 (S27E11 L=309)   368    118        0 (1053/ 458-1157/ 524) -60.3 E12-E08 S30-S23
     10  938  459 (S24E22 L=298)   305    122        1 ( 883/ 435- 991/ 489)   0.0 E26-E17 S25-S23
     11  840  361 (S28E33 L=287)   278     85        1 ( 818/ 321- 856/ 406)  45.0 E35-E29 S31-S25
     12  732  319 (S28E43 L=277)  1349    349        1 ( 598/ 266- 877/ 391) -42.9 E49-E34 S35-S20
  

    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