Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  854 1739 (N58W14 L= 39)   514     92        2 ( 822/1710- 893/1768)  10.3 W09-W20 N54-N59
      2 1497 1585 (N26W51 L= 76)   882    169        1 (1423/1538-1563/1623) -45.0 W45-W54 N21-N29
      3  518 1556 (N54E28 L=357)   391     47        1 ( 486/1536- 541/1581) -15.9 E32-E25 N51-N55
      4 1327 1421 (N22W30 L= 55)   529    151        0 (1261/1379-1376/1441) -39.8 W25-W32 N18-N24
      5 1128 1230 (N16W11 L= 37)   734    163        1 (1064/1170-1183/1275) -76.0 W08-W13 N12-N20
      6  405 1138 (N27E39 L=346)   272     47        1 ( 383/1117- 434/1152)   0.0 E41-E36 N25-N28
      7  832 1086 (N15E10 L= 16)   305     87        1 ( 811/1033- 848/1122)  53.1 E11-E08 N12-N17
      8 1862  905 (S24W58 L= 84)  3520    480        2 (1818/ 794-1928/1061)  54.5 W52-W71 S32-S16
      9  425  997 (N18E39 L=346)   306     70        1 ( 411/ 958- 438/1037)  78.7 E40-E37 N15-N20
     10 1142  931 (S01W04 L= 30)   869    227        3 (1090/ 834-1201/1021)  41.2 W00-W09 S06-N02
     11 1308  686 (S20W09 L= 34)   452    142        0 (1240/ 647-1358/ 708) -40.6 W04-W11 S23-S17
     12 1504  518 (S37W19 L= 45)   685     86        1 (1468/ 471-1535/ 552)  51.3 W16-W22 S41-S35
     13 1213  510 (S28E02 L= 24)   613    182        0 (1138/ 480-1297/ 551) -11.3 E06-W04 S29-S26
     14 1366  393 (S43W06 L= 31)   366     51        1 (1346/ 371-1390/ 421)  11.3 W03-W08 S44-S39
     15 1272  350 (S43E03 L= 23)   349     53        0 (1254/ 325-1294/ 387)  80.5 E03-E01 S45-S39
     16  685  240 (S35E54 L=332)   638     99        0 ( 650/ 206- 735/ 269) -90.0 E56-E51 S39-S32
  

    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