Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  645 1759 (N29E50 L=164)   802    129        3 ( 630/1702- 663/1826) -45.0 E56-E44 N24-N35
      2  496 1613 (N16E52 L=162)   518    111        1 ( 470/1574- 531/1646) -51.3 E54-E48 N13-N18
      3  753 1588 (N20E32 L=182)   394    152        0 ( 705/1535- 796/1642) -24.4 E37-E26 N16-N22
      4  655 1522 (N13E35 L=178)   471    170        1 ( 633/1442- 689/1593) -73.3 E37-E33 N08-N18
      5  382 1489 (N07E55 L=159)   385     88        0 ( 352/1468- 410/1527)  63.4 E56-E52 N05-N09
      6 1753 1445 (N40W44 L=258)   470     64        1 (1730/1412-1778/1476) -78.7 W41-W46 N37-N42
      7  713 1371 (N05E26 L=187)   250     83        2 ( 677/1341- 753/1399) -18.4 E29-E23 N04-N06
      8  470 1169 (S11E37 L=176)   401    133        1 ( 416/1122- 514/1191)  -7.1 E41-E33 S12-S09
      9  165 1067 (S21E66 L=148)   644     75        1 ( 146/1032- 188/1112)  68.2 E68-E62 S23-S18
     10  438  945 (S24E36 L=178)   265     99        0 ( 397/ 935- 495/ 959)  20.6 E39-E31 S25-S22
  

    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