Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1092 1766 (N54E06 L= 86)   588    139        0 (1033/1745-1156/1789)  -4.4 E12-W01 N52-N55
      2  387 1598 (N30E64 L= 28)   600     60        2 ( 368/1570- 403/1636) -45.0 E67-E60 N28-N33
      3 1424 1608 (N45W27 L=118)   400     57        1 (1388/1593-1457/1622)   0.0 W23-W30 N43-N45
      4 1209 1391 (N26W08 L=100)  1893    491        1 (1072/1221-1363/1526) -39.0 E02-W19 N16-N34
      5  308 1452 (N20E62 L= 29)   525     68        1 ( 281/1430- 340/1469)  -6.3 E66-E57 N18-N20
      6 1665 1412 (N32W46 L=138)   528    134        0 (1611/1371-1703/1458) -35.5 W41-W49 N29-N35
      7 1148 1262 (N17W05 L= 97)   351    170        0 (1119/1180-1172/1328) -71.6 W03-W06 N12-N21
      8 1064  708 (S18W06 L= 97)   434    106        2 (1024/ 677-1113/ 745)  38.7 W03-W08 S20-S15
      9 1515  606 (S20W39 L=131)   595    137        2 (1465/ 570-1566/ 627)  33.7 W35-W43 S22-S18
     10  347  547 (S38E53 L= 38)  1210    149        1 ( 287/ 506- 412/ 590)  15.3 E65-E43 S41-S33
     11 1067  190 (S60W23 L=114)   685    121        0 (1044/ 147-1081/ 248) -55.3 W17-W29 S66-S53
  

    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