Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1385 1530 (N34W33 L=355)   830    212        1 (1297/1488-1485/1574) -30.3 W27-W39 N29-N37
      2  797 1538 (N40E13 L=309)   290     64        0 ( 772/1517- 823/1558)  45.0 E15-E11 N38-N41
      3 1490 1338 (N20W35 L=357)   623    181        0 (1448/1255-1541/1401)  36.0 W29-W40 N15-N23
      4 1857 1200 (N07W66 L= 28)   719     61        0 (1827/1177-1884/1225)   5.7 W61-W71 N05-N07
      5 1452 1194 (N12W29 L=351)   249     60        0 (1428/1163-1480/1221) -63.4 W27-W30 N09-N13
      6 1674  856 (S12W43 L=  5)   615    193        0 (1649/ 766-1696/ 951)  76.0 W41-W45 S17-S05
      7  954  499 (S28E09 L=313)   358    127        0 ( 924/ 443- 972/ 545) -90.0 E10-E07 S32-S24
      8  756  437 (S32E25 L=297)   740    195        0 ( 729/ 350- 805/ 511)  47.5 E30-E19 S38-S26
      9  682  295 (S43E39 L=283)   607     78        0 ( 653/ 270- 721/ 325) -59.0 E41-E36 S45-S40
  

    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